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Abstract 

The production of T mesons in fixed target pN collisions is considered. It is shown that 
Non-Relativistic QCD predicts T states to be produced with sizeable transverse polariza- 
tion. The possibility of a measurement of the T polarization at the HERA-i? experiment 
is discussed. 
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1 Introduction 

The Factorization Approach (FA) based on the Non-Relativistic QCD (NRQCD) represents a 
reliable framework to study heavy quarkonium production and decay processes Q. According 
to the FA the inclusive production cross section for a quarkonium state H in the process 

A + B -> H + X (1) 

can be factorized as 

a(A + B^H) = J2^(0\O»\0), (2) 

n m Q 

where the short-distance coefficients, F n , are associated with the production of a heavy quark 
pair in the color and angular momentum state [n]. This part of the cross section involves only 
momenta at least of the order of the heavy quark mass, tuq, and can be calculated perturbatively. 
Two distinct scales are introduced: the heavy quark- ant iquark pair production process occurs 
at small distances, l/(mg), and is factorized from the hadronization phase which takes place at 
large distances, l/(rriQV 2 ). Here v is the average velocity of heavy constituents in the quarkonium, 
with v 2 ~ 0.3 for charmonium and v 2 ~ 0.1 for bottomonium systems. The vacuum matrix 
elements of NRQCD operators, (0|O^|0), describe the evolution of the quark- ant iquark state [n] 
into the final hadronic state H These long distance matrix elements cannot be calculated 
perturbatively, but their relative importance in powers of velocity v can be estimated using 
the NRQCD velocity scaling rules 0. The important feature of this formalism is that the cross 
section of heavy quarkonium production can be organized as double expansion in powers of v and 
a s {mQ). In higher order of v the FA implies that the quark- antiquark color octet intermediate 
states are allowed to contribute to heavy quarkonium production and decay processes. 

Unlike the color singlet long distance matrix elements, each connected with the subsequent 
non-relativistic wave function at the origin, the color octet long distance matrix elements are 
unknown and have to be extracted from the experimental data. The NRQCD factorization 
approach implies universality, i.e. the values of long distance matrix elements extracted from 
different experimental data sets must be the same. However, due to the presently rather large 
theoretical uncertainties §, || and the unknown size of higher twist process contributions || the 
existing experimental data do not allow yet to check the FA universality. 

The data on direct J /if) and if)' production at large transverse momenta at the Tevatron 
indicates that the color octet contribution is dominating. The S state charmonia are produced 
through the gluon fragmentation into the 3 S^' octet state || M. Recent investigations have 
shown that the contribution of color octet states to the charmonium and bottomonium produc- 
tion cross sections is very important at fixed target energies, ^/s ~ 30 — 60 GeV, and reduces 
existing discrepancies between experimental data and predictions of the Color Singlet Model 
(CSM) |S|. However, some experimental data contradict the Color Octet Model (COM) pre- 
dictions. In particular, theoretical predictions disagree with measurements of the polarization 
of J /if) and if)' particles produced at fixed target energies and the COM predicts a too low 
relative yield of the Xd state compared to the Xc2 @ • One possible solution of these discrepancies 
was proposed by Brodsky et al. |J suggesting that higher twist processes, when more than one 
parton from projectile or target participates in the reaction, might give a significant contribution 
to low pr production of J/ if) and x& states. Problems exist also in charmonium photoproduction 
at HERA. The color octet contribution underestimates the inelastic J/ip photoproduction cross 
section at large values of z (z — Ej/^/ E 1 in the laboratory frame) [flCfl . 
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Spin effects in heavy quarkonia production are expected to provide tests for the different 
mechanisms of heavy quarkonium production || f|, [II], |T2| . Predictions for the polarization of 
direct ?/>'s produced at large pt at the Tevatron are free from theoretical uncertainties connected 
with higher twist effects and corresponding measurements will provide an excellent possibility 
to test the NRQCD factorization approach. The observation of opposite sign double spin asym- 
metries in the production of different charmonium states can also be used to discriminate the 
NRQCD FA from the Color Evaporation Model (CEM) which predicts the same spin asymme- 
tries for all charmonium states [JTTJ, jl2| . The CEM also assumes a factorization between the 
production of a heavy quark pair and its hadronization phase. But unlike the NRQCD factoriza- 
tion approach the CEM postulates that multiple soft gluon exchange in the hadronization phase 
destroys the initial polarization of the heavy quark pair and the heavy quarkonium is produced 
unpolarized [13[ . 



At the same time it is not excluded that the mass of the charm quark is not large enough 
to apply the NRQCD factorization approach to charmonium production and decay processes. 
Due to the rather large value of v 2 , about 0.3 for the charmonium system, the Fock states at 
higher order of v 2 may give the essential contribution and can then not be neglected. When 
fitting the values of the long distance parameters the heavy quark spin symmetries are used 
to reduce the number of independent parameters. These relations are valid up to v 2 and may 
get large corrections for charmonium production. In contrast, the NRQCD FA predictions for 
the bottomonium system are more reliable, since the expansion parameter v 2 is much smaller 
(around 0.1), than for the charmonium system. Higher twist processes are also expected to be 
suppressed as A/m b ~ 0.1 (compare with A/m c ~ 0.3). Also the QCD coupling constant is 
smaller for bottomonium system. Therefore, the characteristics of T meson production are more 
appropriate for a correct test of the NRQCD factorization approach. 

In this article we consider the polarization of T mesons produced at fixed target energies. It 
will be shown that in the NRQCD factorization approach T mesons are produced transversely 
polarized, whereas the CEM predicts unpolarized bottomonium production. We present also 
numerical estimates of the projected statistical errors for the measurement of T polarization at 
the HERA-5 experiment at DESY. 

In the next section the bottomonium production in the various subprocesses is discussed. In 
section 3 the T polarization is calculated and numerical estimates for the expected T signal to 
background ratio as well as for the errors of the polarization measurement are considered in the 
section 4. 



2 T Production Subprocesses and Matrix Elements. 

In leading order in a s the different S- and P-wave quark-antiquark states can be produced in 
the following 2^2 and 2^3 subprocesses: 



gluon-gluon fusion 



99- * 8) 99- "Sl^+g 

99-+*p$P gg-^ + g 0) 



gluon-quark scattering 

gq 



3gf > + q gq ^ 3p(i,8) + q (4) 
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quark- ant iquark annihilation 



qq 



35.(8) 



(5) 



where the superscripts (1,8) denote the color singlet and color octet states, respectively. 

The total cross section of T production is given by the sum of direct production cross section 
and the cross section of Xbj states decaying through T mesons: 



o"x = c(T) 



,/ ;r + Br (Xb.j 

J=0,l,2 



TX)a( Xb j) + £ Sr(T(n) - TA>(T(n)). 

n=2,3 



(6) 



The production of each quarkonium state receives contributions from both color octet and color 
singlet states. The relative velocity for the bottomonium system is small, v 2 ~ 0.1, and hence the 
production cross sections can be calculated with only the leading order color octet contribution 
taken into account. 

In table I the color singlet and color octet long distance matrix elements for the production 
of Xbj states are presented. The values of matrix elements are taken from || . The color singlet 
matrix elements are computed from the wave functions of the Buchmiiller-Tye potential tabulated 



m 



3- 



The matrix elements (Of/( 3 Si)) are fitted from the Tevatron data [0]. For the 3P 
bottomonium state the value of this matrix element is obtained by extrapolation from IP and 
2P states §. 



Matrix Elements 


XnO-P) 


X6o(2P) 


Xw(3P) 


W( 3 £i)> 


8.5 ■ 10~ 2 
0.42 


9.9 ■ 10~ 2 
0.32 


0.11 
0.25 



Table I. Matrix elements for the P-wave bottomonia production in units GeV 3 . 

The characteristics of T production at fixed target energies (y/s ~ 40 GeV) differ from 
charmonium production, as it was mentioned in ||. The T(nS) states are mainly produced in 
the quark-antiquark annihilation subprocess through decays of P wave states of bottomonia. 
Unlike J /if) production, the color octet contribution to the T production in the gluon-gluon 
fusion subprocesses is less important. In leading order of perturbative QCD (pQCD) only 
and 3 Pj 8 ' ) color octet states contribute to T production. They are suppressed as v 4 and v 2 
compared to T production through color singlet 3 Si and 3 Pj states, respectively. Hence, due 
to the small value of v 2 ~ 0.1, it turns out that in the gluon-gluon fusion subprocesses the 
main contribution to T production comes from the color singlet states. As can be seen from 
(3), the P-wave bottomonium states are produced in lowest order in pQCD expansion, 0(a 2 ). 
The subprocesses for direct production of S states (T(nS)) are suppressed by the factor a s /n 
compared to those for Xbo and \b2 production. Consequently, the indirect T meson fraction from 
the decay of Xbj states is expected to be large in the gluon-gluon fusion subprocesses. 

The main contribution to P-wave bottomonia production comes from the quark-antiquark 
annihilation subprocess. First of all, the leading color octet and color singlet contributions to 
P-wave quarkonium production scale equally in v 2 , 0(v 5 ), the subleading corrections being only 
of the order 0(v 9 ). In leading order of v 2 , one color octet state contributes to the production of 
Xbj states, namely 3 S[ 8 \ Moreover, in leading order of pQCD only this state is produced in the 
quark-antiquark annihilation subprocesses. The values of the matrix elements (Og bJ ( 3 Si)) are 
larger as compared to the color singlet matrix elements for corresponding P states, connected to 
the quarkonium wave functions derivatives at the origin (see table I) 0. Furthermore, at fixed 
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target energies the qq luminosity effectively increases compared to the gg luminosity due to the 
large mass of the (bb) system. 

The T production cross sections in gluon-gluon fusion and quark- ant iquark annihilation sub- 
processes are presented in table II. The 6-quark mass of = 4.9 GeV is chosen as in ref. |7j] in 
order to extract the values of color octet long distance parameters from the Tevatron data for the 
T production. The cross sections are calculated using the GRV LO parton distribution functions 
|i~5|| evaluated at the factorization scale Q 2 = 4m\. The long distance color octet parameters 
for direct T(nS) production are taken from ||. For the decays of Xbj{3n) states[| to T(3S) the 
same branching ratios are assumed as for the corresponding n = 2 states |J . 



subprocess 


gluon-gluon fusion 


quark- ant iquark annihilation 
through Xbj 


through T(nS) 


through Xbj 


<r(T) 


0.013 nb 


0.016 nb 


0.24 nb 



Table II. T(IS') production cross sections in different subprocesses. 

From table II one finds that the direct T production cross section for the adopted values of 
the long distance parameters is more than one order of magnitude smaller than the total cross 
section. Even in the color singlet model the T mesons are mainly produced through decays of Xbj 
states. On the other hand, the color octet contribution seems to be dominant in Xbj production. 
The main contribution comes from the 3 Si state, produced in quark- ant iquark annihilation. 

As was already shown in ||, the color octet contribution reduces the large discrepancy 
between the CSM prediction for the total T production cross section and experimental data 
16, O, 111]. Nevertheless, there remain large uncertainties due to contradicting experimental 



results. The cross sections obtained by integration of xp distributions for T(1S) production 
presented in [|l6j and |17] are three and four times smaller, respectively, than the central value 
quoted in JT£[] , 270 pb/nucleon. In addition, the theoretical value of the cross section strongly 
depends on the assumed mass of the 6-quark. It is therefore impossible to extract the color octet 
matrix elements with reasonable accuracy from the fixed target bottomonia production data. 



3 T Polarization 

As it was already mentioned in , the measurement of the cross section for direct and indirect 
production of T's would provide crucial information about the color octet mechanism, e.g., new 
constraints on the long distance color octet parameters would emerge. Such a measurement 
requires the reconstruction of Xbj states in the T + 7 decay mode which is not a trivial task at 
fixed target experiments due to the small transverse momentum of the emitted photon. 

Another possibility to check the NRQCD factorization approach is to measure the T po- 
larization. In H — > £ + £~ decays the polarization of 5-state quarkonium is determined by the 
polar-angle distribution of its decay leptons with respect to the beam direction in the meson rest 
frame. Integrating over the azimuthal angle the distribution has the form 

oc 1 + a cos 9, (7) 



dcosO 



2 As was suggested in ||, the large yield of the T(3S f ) state at the E772 experiment at FNAL (T(| could be 
explained by assuming the unobserved Xf>j(3P) states to lie below the open beauty threshold. 
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where 9 is the angle between the positively charged lepton, £ + (£ — e, //), and the beam axis in 
the quarkonium rest frame. The parameter a in the angular distribution can be related to £, 
the fraction of longitudinally polarized T mesons: 

1 - 3£ f 1 for £ = 



1+^ l-lforC = l KJ 

The calculation of J/^ and if)' polarization at fixed target energies and Tevatron collider 
energies was performed in j|, |9|, [L^, |2(J . The most general method to calculate the cross sections 
for heavy quarkonium production with definite polarization within the NRQCD factorization 
approach was proposed by Braaten and Chen ||21|| . Tang and Vanttinen used the covariant 



projection method to calculate cross sections for polarized J /if) and if)' production p0| . The 
XbjinP) states are produced in quark-antiquark annihilation through only one color octet state, 
3 S[ 8 \ For this particular case both methods give the same result. We used the formulae derived 
in [|(J to obtain the polarization of T(nS) states produced in cascade, qq — > 3 S[ 8 ^ — > Xbj — > 
T(nS)+ r - 



a(qq - bb( 3 S[ 8) ) - Xbl + 9 - T + 7 + g) 

'%\ - M 2 /s)Br( Xlb -> T + 7 )(0 8 ^(3 5l ))^_^£, ( 9 ) 



27M 5 

°{qq - bb( 3 S{ 8) ) ^ Xb2 + g ^T + 1 + g) 

= - MVs)Br( Xb2 — > T + T)(0f 2 ( 3 ^)) 47 ~ 2 2 o HA ° - (10) 

The scalar x^o state yields unpolarized bottomonium S"-wave states. It follows from (9) and 
(10) that Xbi an d X b2 intermediate states lead to values of a = 0.2 and a = 0.29, respectively. 
Taking into account all transitions from X bj(lP) and X bj(2P) states to T(1S) we obtain for the 
polarization parameter in quark-antiquark annihilation subprocess 

a ~ 0.24. (11) 

This result remains practically unchanged if we add contribution from X bj{3P) states with the 
same branching ratios as for the corresponding 2P states. It is worth mentioning that the value 
a ~ 0.24 represents only a lower theoretical limit for the polarization. This value is calculated 
taking into account only the quark-antiquark annihilation subprocesses. In gluon-gluon fusion 
subprocesses the polarization of T mesons is larger due to the dominant contribution from 
X b2 decays which yield pure transverse polarization p. For the values of octet long distance 
parameters (table I) the size of the polarization is a ~ 0.3 and hence also exceeds 0.24. However, 
large uncertainties in color octet long distance parameters for P-wave bottomonia does not allow 
to compute the relative importance of the various subprocesses, so that the lower bound of 
a = 0.24 remains as the only firm prediction. 

In contrast to charmonium production the higher twist effects for the bottomonium system 
are expected to be small. In particular, the higher twist effect suggested in 0, when more than 
one parton from the projectile or target is involved in the heavy quarkonium states production, is 
expected to be negligible for T production. Two partons should be within a transverse distance 
of 0(l/mg) in order to interact with the other parton and to produce a heavy quark-antiquark 
bound state. Consequently, such a higher twist process is suppressed by a factor of 0{K 2 Q CD /mQ) 
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H . To explain the discrepancies between the CSM predictions and measured relative production 
rate Xci/Xc2 it is suggested that the above suppression can be compensated by a kinematical 
enhancement ||. Thus the Xci production cross section in the higher twist process is expected 
to be at the same level as the Xc2 production cross section from gluon-gluon fusion subprocesses. 
In bottomonium production such a mechanism is suppressed by the mass of the bottom quark 
0(Ag CD /m^) , i.e. the corresponding cross section should be one order of magnitude smaller than 
in the charmonium case. The effect of the kinematical enhancement is also smaller due to the 
larger mass of the bottomonium system. Moreover, as can be seen from table II, the gluon-gluon 
fusion subprocesses give small contributions to T production at fixed target energies compared 
to J/ip production. Therefore, the measurement of T polarization at fixed target energies will 
allow to distinguish between the NRQCD approach and the CEM, which predicts unpolarized 
production of all bottomonium states. 



4 Expectations for HERA-B 

HERA- 1? is an experiment presently set up at DESY which uses the HERA 920 GeV/c proton 



beam incident on various nuclear targets |P2jj . Being optimized for the study of the various 
aspects of .B-physics, i.g. the measurement of the CP asymmetry in the — > J/tpK® decays, the 
experiment is well suited to perform accurate T polarization measurements due to the following 
reasons: 

• large acceptance of the apparatus (polar angle coverage in the laboratory frame from 
10 mrad up to 250 mrad), 

• good momentum/mass resolution which allows to separate the T(nS) states, 

• high statistics due to high interaction rates (40 MHz). 

In the current analysis the acceptance for T production and the mass resolution in the T mass 
region are determined using a detailed HERA-5 detector simulation P5[ . The T mesons are 
generated with the PYTHIA 5.7 event generator |24j] and then processed through the simulation 
of the full detector. The cross sections for T production are taken from the measurements of the 
E605 experiment (^/s ^ 38.8 GeV) |] 



Br-^\ y=0 (T(lS) + T(2S) + T(35)) = 2.31 pb/nucleon; 
ay 

Br^(T(2S)) Br^(T(3S)) 

— ^L— — — — = n ^1 ■ dy — — n no 

Br%{T{lS)) U - di ' Br%(T(lS)) Um ^ 

As it was noted above (in section 2), this experiment gives the lowest lying value for the T 
production cross section , about 130 pb/nucleon for cr(Y(lS') + T(2S') + T(3S)). 

In the following only the muonic decay channel is considered because the expected mass 
resolution is better. Figure 1 shows the T — > mass distributions on top of the Drell-Yan 

pair production expected in one year (10 7 s) of HERA-I? running with 40 MHz interactions rate. 
Cut on the muon momenta, p > 5 GeV/c and pt > 0.5 GeV/c, are applied to satisfy the trigger 
requirements. The expected statistics are about 10000 fully reconstructed T events per year. 
The background is largely dominated by the Drell-Yan process. It is simulated by PYTHIA with 
the cross section calculated at the leading order, thus a K factor of 2.3 is taken into account 



according to |T7|]. As can be seen from Fig.l the predicted mass resolution of about 1% and a 
signal to background ratio S/B~ 1.3 (in a ±2<r mass window around the T(1S*)) allows to clearly 
observe the T mass peaks above the background and to separate the T(IS') and T(2S) states. 
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Figure 1: T — * fi + fi~ mass distribution in 
one year running of HERA- 1? 
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Figure 2: Acceptance for the T events as 



a function of the CMS polar angle 



Figure 2 shows the acceptance for T events as a function of the cosine of the polar angle 
between the positive muon momentum and the beam direction in the rest frame of T. In almost 
the full range of cos 9 the acceptance is close to unity. Large acceptance corrections are expected 
only for | cos#| > 0.8. 

The accuracy of the polarization measurement is estimated by simulating various initial 
polarizations for T mesons and Drell-Yan pairs with various signal to background ratios. After 
corrections for the acceptance the resulting angular distributions are fitted with the function: 

dN 

- = S{l + acos 2 6) + B(l + a'cos 2 6). (13) 

a cost' 

Here S and B are signal and background normalization factors and a and a' determine corre- 
sponding polarization parameters. This procedure allows to extract the polarization not only 
of T mesons but also of the Drell-Yan pair, which is an interesting issue by itself. Assuming 
that Drell-Yan polarization and signal to background ratio could be precisely measured in the 
HERA-5 experiment, the number of free parameters in the fitting function is reduced to 2. 
Results on the parameter a obtained from the fit for the different MC input values and for 10000 
reconstructed T events are shown in table III. 



Input values 
of a 


S/B=1.3 


S/B=l 


S/B=2 


a' = 


a' = 1 


a' = 1 


a = 0.25 


0.21 ±0.08 


0.21 ±0.08 


0.22 ±0.09 


0.23 ±0.08 


a = 0.5 


0.52 ±0.10 


0.48 ±0.09 


0.51 ±0.09 


0.54 ±0.09 


a — 1 


0.91 ±0.12 


0.96 ±0.11 


0.96 ±0.12 


0.98 ±0.10 



Table III. Parameter a obtained from the fit for different MC input values. 

We note that the expected statistical error on a is largely dominated by the T statistics 
rather than by the value of the background polarization and the signal to background ratio. 
Figure 3 illustrates the statistical error 5a as a function of the number of reconstructed T events 
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for various cv's. As can be seen, an accuracy on the 
polarization parameter 5a ~ 0.08 for a = 0.25 can 
be achieved for one year of the HERA-5 running. 
The expected error will allow to distinguish between 
a = (CEM) and a ~ 0.24-0.3 (NRQCD FA) with 
a 3cr significance. 

5 Conclusions 

The polarization of T mesons is calculated at fixed 
target energies (y/s ~ 40 GeV). It is shown that in 
the NRQCD factorization approach T mesons are ex- 
pected to be produced transversely polarized; the pa- 
rameter a for the polar angle distribution of quarko- 
nium decay products is about 0.24 -=-0.3. In contrast, 
the color evaporation model postulates that multiple 
soft gluon exchange in the hadronization phase de- 
stroys the initial polarization of the heavy quark pair 
and quarkonium is produced unpolarized. 

Higher twist effects are expected to be small due 
to the large mass of the 6-quark. The contribution of higher Fock states in bottomonium produc- 
tion are more suppressed than in the charmonium case, the relative velocity for the bottomonium 
family is about v 2 ~ 0.1. Thus the measurement of the T polarization provides an excellent 
opportunity to test different mechanisms of heavy quarkonium production. In particular, it 
allows to distinguish between the NRQCD FA Q and the CEM [Ej|. On the other side, the 
observation of an extremely large polarization will indicate that T mesons are mainly produced 
through color singlet states and that the color octet parameters for Xbj production extracted 
from the Tevatron data should be much smaller than known at present. 

The Monte Carlo simulation shows that the projected statistical error for the measurement 
of the polarization parameter a is about 0.08 for a ~ 0.25 in one year running of the HERA- 
B experiment. The simulation is done only for the decay channel. A statistics gain of 

almost a factor 2 is expected from the e + e~ channel. The simulation conservatively assume the 
lowest value for the T production cross section measured in different experiments [in| [17], [Tj| 
at y/s ~ 38.8 GeV which corresponds to a HERA proton beam momentum of 800 GeV/c (the 
current value for HERA is 920 GeV/c). 
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